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THE  EXPANSION  OF  A  RAREFIED  GAS  INTO  A  VACUUM 

P.  Molmud 

Space  Technology  Laboratories,  Inc.,  Los  Angeles  45,  California 

ABSTRACT 

The  expansion  of  an  initially  uniform  gas  sphere  into  a  vacuum  is 
described.  The  analysis  px-oceeds  on  the  basis  of  free  molecular  flow. 

It  is  demonstrated  that  collisions  will  not  seriously  modify  the  expan¬ 
sion  process  as  long  as  the  initial  mean  free  path  of  a  molecule  within 
the  sphere  is  much  greater  than  the  initial  diameter  of  the  sphere. 

Expressions  are  derived  for  the  momentum  and  energy  flow  as¬ 
sociated  with  the  expansion  of  such  a  sphere. 

The  expansion  of  a  gas  which  fills  half  space  and  cdso  of  a  gas 
initially  of  the  form  of  a  uniform  cylinder  are  described. 

The  interaction  of  an  expanding  gas  cloud  with  a  specularly  re¬ 
flecting  surface  is  determined. 

The  interaction  of  an  expanding  gas  cloud  with  a  diffusely  reflecting 
surface  is  considered. 

I  INTRODUCTION 

Let  us  consider  a  rocket  vehicle  operating  in  interplanetary  space.  The 
vehicle’s  rocket  motor  may  be  operating  under  steady-state  thrust  or  sporadi¬ 
cally  in  short  bursts,  either  for  orientation  or  for  braking  purposes.  We  are 
concerned  with  the  phenomena  that  may  be  associated  with  the  release  of  exhaust 
gases  from  the  vehicle.  For  example:  the  released  gases  may  by  their  very 
expansion  into  the  vacuum  react  on  the  vehicle.  The  release  gases  may  also, 
iunder  the  influence  of  solar  radiation,  fluoresce  and  ionize.  Furthermore 
some  of  the  ions  so  produced,  or  already  present  m  the  exhaust  gases,  will  be 
trapped  by  the  local  magnetic  fields. 

In  order  to  predict  the  magnitude  and  spatial  extent  of  such  effects,  it  is 
necessary,  at  the  very  least,  to  be  able  to  predict  the  temporal  evolution  of 
the  released  gases. 
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In  this  paper  we  consider  the  expansion  of  a  short  burst  of  gas  (a  puff) 
into  a  vacuum.  The  continuum  treatment  of  this  apparently  simple  problem 
IS  amazingly  difficult  and  by  no  means  complete.  ^ 


^See,  for  example,  the  treatment  by  J.  B.  Keller,  Quarterly  of  Applied 
Mathematics,  14,  (1956). 


Therefore  we  confine  ourselves  to  a  treatment  of  this  problem  which  can  be 
carried  through  to  completion  and  where  the  region  of  validity  of  the  solution 
can  be  fairly  well  defined.  The  treatment  we  employ  is  that  of  the  free 
molecular  flow  of  gases;  that  is,  we  assume  that  the  molecules  in  this  puff 
of  gas  move  to  infinity  without  colliding  with  one  another. 

The  solution  is  simple  and  proceeds  in  a  straightforward  manner.  After 
we  obtain  the  solution,  we  investigate  the  consequences  of  assuming  a  finite 
collision  cross-section  for  the  molecules.  We  find  that  the  assumption  of 
free  molecular  flow  is  not  inconsistent  with  the  assumption  of  the  finite 
collision  cross-section,  provided  the  in...ial  diameter  of  the  puff  is  much  less 
than  the  initial  mean  free  path  of  the  molecules. 

II.  THE  FREE  MOLECULAR  EXPANSION  OF  A  PUFF 
OF  GAS  INTO  A  VACUUM 

We  shall  simulate  our  puff  of  gas  by  a  spherical  bubble  of  gas  in  containment 
by  some  sort  of  skin  which  can  be  made  completely  permeable  at  will.  That  is, 
the  skin  will  retain  the  gas  until  a  signal  is  given,  at  which  point  the  skin  ceases 
to  exist  and  the  gas  is  allowed  to  expand. 

We  shall  furthermore  assume  that  the  gas  molecules  initially  have  a 
Maxwellian  distribution  of  velocities  at  a  temperature  T  ,  and  do  not  collide  with 
each  other  during  the  expansion  (zero  collision  cross-section). 
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1.  The  Source  Function  for  Free  Molecular  Expansion  into  a  Vacuum 


We  asstime  that  the  skin  bounding  the  sphere  of  gas  disappears  at  t  =  0. 
Since  the  particles  do  not  interact  with  each  other,  it  is  sufficient  for  the 
complete  description  of  the  problem  to  know  the  initial  vector  velocity  of  each 
particle.  Then  we  can  allow  these  particles  to  travel  in  straight  lines  after 


t  =  0. 

Let  us  take  an  element  of  volume  d  t  containing  d  N  particles  in  the 
original  sphere.  Since,  by  hypothesis,  these  particles  nave  a  Maxwellian  velocity 
distribution,  we  may  immediately  write  the  fraction  f,  of  the  dN  particles 
with  radial  velocities  ranging  between  v  and  v  +  dv.  Thus: 

^  v^dv  (1) 

where  p  =  m/2KT  (m  is  the  mass  of  the  molecule,  K  is  the  Boltzman  con¬ 
stant  and  T  is  the  temperature  of  the  gas). 

At  a  time  t,  after  the  start  of  the  expansion  these  particles  with  velocities 
V  to  V  +  dv  are  contained  between  two  concentric  spheres  having  a  center  at 
dr  and  radii  r  =  vt  and  r  +  d  r  =  (v +  d v)  t.  Therefore,  the  density  of  these 
particles,  dp,  is  given  by 


,  fdN 

dp  =  - . 

4Trdvt  V 


(2) 


Combining  (2)  with  (1)  and  setting  v  =  r/t,  we  obtain 


3/2  ,-P(rW) 


dN 


(3) 


Equation  (3)  presents  the  density  everywhere  in  time  and  space,  due  to  the 
free  expansion  of  a  point  Maxwellian  ensemble  of  dN  particles.  It,  therefore, 
deserves  the  title  of  "source  function  for  free  molecule  expansion." 
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This  source  function  is  reminiscent  of  the  source  functions  for  diffusion 


and  thermal  conductivity.  There  is  an  important  difference,  however,  and  that 
is,  time  is  to  the  first  power  in  the  exponent  and  to  the  3/2  power  in  the  denomi¬ 
nator  of  these  latter  functions. 

As  can  readily  be  demonstrated. 


flOO 


4  ITT  dpdr  =  dN 


independent  of  time;  and,  also,  the  time  integral  of  the  flux  of  particles  through 
the  surface  of  a  surroundmg  sphere 

100 


f  .  2  ^  d  dt 

J4nr  —dpdt  , 

o 


is  equal  to  dN,  independent  of  r.  These  are  necessary  properties  of  a  source 
function  and  are  possessed  by  (3). 

2.  The  Source  Function  for  a  Moving  Source 

The  source  function  (3)  can  be  generalized  to  take  into  account  the  motion 
of  the  source  (i.  e.  ,  the  initial  Maxwellian  distribution  is  centered  about  a  drift 
velocity  V). 

Let  us  treat  the  case  of  the  moving  observer,  since  the  moving  source  can 
always  be  transformed  to  such  a  case.  We  have,  according  to  Figure  1,  an 
observer  moving  with  velocity  V.  At  t  =  0  he  is  a  distance  r^,  from  the  source, 
8,  and  V  makes  an  angle  of  0  with  respect  to  the  radius  vector,  r^,  drawn 
from  the  source. 
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Now,  the  radial  distance  from  the  source  will  be  given  by 


r  =  i/r^  +  V^t^  -  2Vtr  cos  0 
Vo  o 


Thus  since  solution  (3)  is  valid  for  all  r  and  all  time,  the  density 
encountered  by  the  moving  observer  is: 


=  11 


-p(r^  +  V^t^-2Vtr^  cos  e)/t^ 


The  situation  as  viewed  by  a  stationary  observer  and  moving  source  is 


thus: 


-p(r^  +  V^t^  +  2Vtr^  cos  e)/t^ 


dp  = 


When  the  source  is  in  motion  V  can,  in  general,  be  a  function  of  position  so 
that  (5a)  is  the  source  function  for  a  body  of  gas,  the  individual  parts  of  which 
are  in  relative  motion  (in  the  hydrodynamic  sense)  at  t  =  0. 

3.  The  Uniform  Gas  Sphere  in  Free  Expansion 


/A 


Figure  2 


Using  the  source  (3)  we  may  now  set  up  an  expression,  valid  everywhere 
in  space  and  time,  for  the  density  produced  by  the  free  expansion  of  an  initially 
uniform  sphere  of  gas  with  initial  density  p^. 
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With  reference  to  Figure  2; 


,a  ^Tf 


p(R,t)  = 


O  ITT 


2  TT r  sin  0d  6d  r 


(6) 


o  u  o 


2  2  2 

where  =  dN/dr.  Since  i  =  r  +R  -  2rR  cos  6  we  may  peiform  the 
tedious  integration  and  obtain 

p(R,t)  =  ^|Erf  j{R  +  a)-^j- Erf  [(R-a)^j| 

2./.2 


+  [e-  (a  +  R)^  p/r  _  (a  -  R)^  p/t^ 

2R^L  J 


(7) 


where 


Erfx 


dy  . 


Equation  (7)  is  simplified  at  either  R  s  0,  the  center  of  the  sphere,  or  R  >>  a, 
where  a  is  the  initial  extent  of  the  sphere. 

At  R  ~  0  our  solution  becomes 


p(O.t)  = 


(8) 


pa  fliT  2 

and  for  large  R,  (R>>  a),  i«R  and  J  J  Pn  ^irr  sin  0d0  =  N,  the  total 

J o  Jo 


number  of  particlss  in  the  original  sphere.  Thus  for  large  R 


p(R,t)  =  I'  ® 


N 


(9) 


p  having  exactly  the  same  form  as  the  source  function. 

We  can  cast  Eq.  (9)  into  a  form  with  more  physical  meaning  by  borrowing 

2 

a  description  employed  by  J.  B.  Keller.  We  assume  the  gas  to  have  the  specific 


2 

J.  B.  Keller,  "On  the  Solution  of  the  Boltzmann  Equation  for  Rarefied  Gases." 
Comm.  Pure  Appl.  Math.,  1,  275  (1948). 
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heat  ratio  y  and  the  velocity  of  sound  in  the  original  sphere  is  = 
Then  p  =  ylZC^  and  therefore: 


fy  (KT/m). 


-  i  »■ 


(C„t)- 


An  inspection  of  Eq.  (9a)  readily  reveals  that  the  main  body  of  the  gas  is 
contained  in  a  sphere,  the  radius  of  which  is  expanding  with  the  velocity  . 

4.  The  Pressure  and  Impulse  Due  to  an  Expanding  Gas  Sphere 

The  pressure  at  any  point  uu :  to  an  expanding  gas  sphere  may  be  determined 

by  computing  the  momentum  transport  through  a  unit  area  per  imit  time.  We  set 

up  a  unit  area  as  in  Figure  Z,  perpendicular  to  the  radius  vector,  R.  The 

normal  momentum  flow  per  unit  area  per  unit  time  due  to  each  source  element  is 
2  2 

(dp)mv  cos  Therefore,  the  pressure  at  R  is: 


P(R.t)  = 


QUO 


T  j  m  v^cos^i|<  r^sin0d0dr 


.  (10) 


2  2  2 

where  v  =  i  /t,  i  cos  =  R  -  r  cos  0,  i  sr  +R  -  2rR  cos  0. 

2 

Again,  at  large  distances  (R^>  a),  i  ^  R,  cos  (j;  and  P  becomes 


M  r"  I' 

1  — 5  ® 

r  {C^tr 


(lOa) 


where  M  =  Nm,  the  total  mass  within  the  original  sphere. 

The  momentum  or  impulse  transmitted  through  the  unit  area  is  given  by 
the  time  integral  of  (10).  The  total  impulse  per  unit  area  is: 
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At  large  distcinces,  this  reduces  to 


I 


MC 

o 


(lla) 


This  last  result  is  readily  obtained  by  other  means:  Consider  dN  particles 
in  a  Maxwellian  distribution.  The  momentum  contained  in  those  particles  which 
have  velocities  v  to  v  +  d  v,  in  a  solid  angle  dl2  is 


dl  =  m  v' 


(?) 


e  ^  dvdfidN  . 


(12) 


The  average  momentum  in  the  direction  of  the  infinitesimal  solid  angle 
d  f2  is  obtained  by  integrating  over  the  velocity,  thus 

.3/2 


dl 


=  (!)■ 


m  d  N 


df2 


(12a) 


Now  instead  of  dfl;  we  can  substitute  its  measure,  the  area,  d  A,  subtended 
by  this  solid  angle  at  the  distance  R,  divided  by  Th..s 


d  1  m  d  N 


(12b) 


equivalent  to  (lla). 

5.  Total  Power  and  Energy  Transmitted  through  Unit  Area  by  an  Expanding 
Gas  Sphere 

The  flux  of  energy  is  given  by  1/2  m  v^  (dp)  v  cos  ijj  (refer  to  Fig.  2)  where 
V  =  i/t.  Thus,  the  power  transmitted  through  a  unit  area  is  given  by  the 
following  expression: 


W 


2  2 

e  ^  cos  sin  0  d  0  d  r. 


(13) 


Again  at  large  distances  (R  p->a),  W  reduces  to 
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(13a) 


W 


M 


r3  e- 


where  M  =  Nm,  the  total  mass  in  the  initial  sphere. 

The  total  energy,  £,  passing  through  a  unit  area  is  obtained  by  integrating 
£q.  (13a}  with  respect  to  time.  The  result  is 


£ 


M 

~~2 — 
R  TTp 


(14) 


Thus  the  total  energy  crossing  a  large  sphere  about  the  original  one  is  just 
2  3  M 

4irR  times  £q.  (14)  ^  •  This  is  a  resonable  result,  for  the  average 

3  m  3 

energy  of  the  molecules  originally  in  the  sphere  was  ^  *p~  or  KT.  Thus 

3  M 

the  total  energy  in  the  sphere  was  ^  ~  and  this  total  energy  is  transported 
through  the  boundaries  of  the  large  bounding  sphere,  yielding  £q.  (14). 


III.  TH£  CONSIST£NCY  CRIT£RIA  FOR  TH£  VALIDITY 
OF  FRE£  MOLECULAR  FLOW 

We  have  up  to  this  time  assumed  that  the  molecules  do  not  interact  with 
one  another,  that  is,  they  have  zero  collision  cross-section.  We  have  des¬ 
cribed  the  expansion  of  a  sphere  of  gas  imder  these  assumptions.  We  now 
ask  the  question:  Under  what  criteria  is  this  description  vedid  for  real 
molecules  with  finite  collision  cross-sections  ? 

We  proceed  to  determine  these  criteria  by  the  following  method:  Each 
molecule  in  the  gas  is  to  be  endowed  with  a  finite  collision  cross-section, 
but  each  molecule  is  artificially  constrained  to  move  in  a  straight  line  to 
infinity,  even  though  it  might  collide  with  another  molecule.  The  expansion 
of  the  gas  is  described  according  to  the  methods  of  Section  II,  but  now  we 
determine  the  number  of  collisions  the  average  molecule  encounters  on  its 
rigid  and  prescribed  path  to  infinity.  If  this  number  is  much  less  than  one, 
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then  we  are  correct  in  assuming  that  most  of  the  molecules  escape  the  sphere 
without  suffering  a  collision.  Consequently,  the  description  of  the  expansion 
according  to  free  molecular  flow  is  valid. 

1 .  The  Collisions  Between  one  Particle  from  one  Source  and  a  Swarm  of 
Particles  from  Another  Source 

We  first  consider  two  point,  particle  sources,  Sj  and  s^,  and  a  particle 
of  velocity  V,  emanating  from  s^  at  t  =  0  and  proceeding  along  a  straight  line. 
We  shall  assume  that  this  particle  does  not  collide  with  any  other  s^  particles. 
The  Sj  particle  will  pass  through  a  swarm  of  particles  emanating  zilso  at  time 
t  =  0  from  &2  (see  Fig.  3). 


Figure  3 

The  collision  frequency  between  the  Sj  particle  and  the  s^  particles  is 
given  by  P2®^12’  ^2  density  of  the  82  particles  encountered 

by  the  first  particle  in  its  path,  Q  is  the  collision  cross-section  and  Vj2  is 
the  relative  velocity  between  the  two  particle  species.  When  the  expression 
P2QV12  is  integrated  over  all  time,  the  result  will  be  the  total  number  of 
collisions  suffered  by  the  Sj  particle  in  traveling  an  infinite  distance  from 
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The  relative  velocity,  between  particle  and  the  particles  it 

meets  on  its  rectilinear  path  is  just  r/t,  independent  of  the  velocity  of  the 
Sj  particle.  This  is  made  evident  by  first  considering  those  s^  particles  of 
zero  velocity.  The  s^  particles  they  meet  have  a  velocity  r/t.  Now  consider 
Sj  particles  with  velocity  V  in  any  direction.  The  s^  particles  that  are  going 
to  be  encountered  must  have  this  velocity,  V,  vectorially  added  to  r/t.  Thus, 
the  relative  velocity  remains  r/t. 

Let  us  denote  dn  as  the  total  number  of  collisions  encountered  by  a 
^.Article  of  velocity  V  in  the  presence  of  source  s^.  Thus: 


dn 


(15) 


but 


'  ‘’o  (1)  ^ 


3/2 


2  2  2  2 

where  dr  is  the  volume  of  the  source  3^,  its  density  and  f  =r  +V  t  - 


2rVtcos  0,  l^cf.  Eq.  (5)J.  Thus: 


dn  =  Qrp_jdT^-|j  J 


pV^zcose 


(15a) 


2.  The  Average  Number  of  Collisions  between  Particles  from  one  Source 
and  a  Swarm  of  Particles  from  Another  Source 

Now  the  number  of  particles  from  s.  with  a  velocity  V  to  V  +  d  V,  at 

^  2 

an  angle  0  and  in  a  solid  angle  Zv  sin  0  d  0,  is  2tt  (P/tt)^^^  sin  0  d  0  e’^  ^  d  V. 

We  now  weight  (15a)  with  this  factor  and  integrate  over  all  angles  and  velocities, 
thus  obtaining  the  average  number  of  collisions  encountered  by  all  particles 
emanating  from  Sj  and  proceeding  to  infinity. 
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We  perform  the  0  integration  first  and  obtain 


<dn>.  = 


e  ^ 


-pv‘ 

V  ze 


Next,  the  velocity  integration 


00 


<^dn>  2 

Zr  ir 


=  MN  1  /  ,azr_l 


_{z  -  1 )  +  1  (z  +  1  ) 


-^1 

+ 1 )  +  1 J 


(17) 


where  dN  =  P^dr  and,  finally  the  z  integration 


<dn  > 


e,  V 


QdN 

1 

Zttt 


(18) 


3.  The  Average  Number  of  Collisions  Encoimted  by  a  Particle  in  Escaping 
From  a  Uniform  Spherical  Gas  Cloud 

We  now  determine  the  average  number  of  collisions  encountered  by 

particles  from  s^  when  s^  is  located  in  an  initially  uniform  spherical  cloud 

and  8,  is  a  distance  R  from  the  center. 

1  o 

We  choose  a  polar  coordinate  system  such  that  s^  is  on  the  polar  axis, 
as  in  Figure  4. 


Figure  4 
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2  2  2 

Then  r  =  R  +  R  -  2RR  cos  8  and  therefore  from  (18) 

o  o 


a  ? 


<n> 


e.  V 


=  Qp, 


‘'o  “o 


R  dR  sin  8  d9 

R^  +  R^  -  2RR  cos  8 
o  o 


(19) 


Performing  :he  8  integration  first,  we  obtain 


< 


Qp  ^  fR  +  K  \  ^ 

«>8.v  -  nr  J  • 


(19a) 


Then  the  R  integration: 


Qp, 


<n> 


8,  V  2 
Thus,  when  R  =  0,  <n> 


2  /a+R 

ir  '"IrTir 

o  V  o 


<a-  R. 
"*■  ^"0  a  +  R 


+  2  a 


(19b) 


8  V  =  2QPo^‘ 


e,v=  “'’o*  “■* 


when  e  ^  Qp^a. 

Since  P^Q  1^  the  reciprocal  of  the  mean  free  path  and  2  a  is  the  diameter 

of  the  sphere,  we  see  that  the  requirements  for  neglecting  collisions  are  pre> 

cisely  what  we  would  expect,  namely,  that  the  diameter  of  the  sphere  be  much 

less  than  the  mean  free  path  of  the  moleculefi,  or  equivalently  that  . 

•  1 6  2 

Since  the  cross  section  is  usually  on  the  order  of  10’  cm  this  means  that 
1 6 

p^<^<  10  /2a.  Thus,  for  example,  the  expansion  of  a  sphere  initially  with  a 

1 6  +3 

diameter  of  one  millimeter  and  with  a  particle  density  of  10  per  cm  may 


be  treated  as  in  the  first  two  sections  of  this  report. 
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APPENDIX  I 


OTHER  TRANSIENT  PROBLEMS  INVOLVING  THE  FREE 
EXPANSION  OF  G.4SES  INTO  A  VACUUM 

Using  the  source  function.  Equation  (3),  we  may  solve  other  problems 
involving  the  free  expansion  of  gases  into  a  vacuum. 

1 .  The  Free  Expansion  of  a  Gas,  Filling  Half  Space,  into  a  Vacuum 

We  consider  the  problem,  already  solved  by  Keller^  of  the  expansion  of 
a  gas  for  which  p  =  for  x<  0  and  p  =  0  for  x^O  at  time  t  4  0. 

3 

cf.  footnote  2. 


First  we  determine  the  source  fvmction  of  a  flat  slab  extending  to  infinity 
in  the  ±y  and  ±z  directions  but  of  thickness  dx.  Consider  Figure  I>1  where 


dy 

y 


dx 

Figure  I-l 

we  have  the  infinite  slab'of  density  p^  and  an  observer  at  distance  R  -  ' 

normal  to  the  slab.  The  element  of  volume,  dr,  is  Z^rydydx  and  also 
2  2  2 

i  =  (R-x)  +y  .  Now  using  the  source  function  (3),  we  find  the  density 
at  R  due  to  the  initially  uniform  infinite  slab. 


dp 


y  d  y  d  X 


(I.  i) 
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The  source  function  (1.  1)  when  integrated  over  x  from  -oo  to  0  will  give  us 
the  solution  to  the  problem.  Thus: 


p(R.T) 


^-p(R-x)^/t^ 


-00 


dx  . 


Let  p/t^  (R  -  x)^ 


and  we  obtain 


(1.2) 


(I.  3) 


the  same  result  obtained  by  Keller.  This  solution,  as  Keller  points  out,  also 
describes  the  expansion  of  a  semi-infinite  gas  in  an  infinite  cylinder  of  rectarngular 
cross-section  with  specularly  reflecting  walls. 


2.  The  Free  Expansion  of  a  line  Source  of  Gas 

Consider  a  line  source  of  gas  of  infinite  length  and  of  linear  gas  density, 
dN/dl  .  Then  consider.  Figure  1-2,  which  illustrates  the  computation  of  the 
density  at  time  t  at  a  distance  R  from  this  line. 


Now  r^  =  f  ^  +  R^  and  from  Eq.  (3) 


7 


(1.4) 


We  use  this  result  immediately  in  the  next  section. 
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3.  The  Free  Expansion  of  a  imiform  cylinder  of  Gas,  of  Infinite  Length, 
into  a  Vacttum 

We  determine  the  source  function  of  a  cylindrical  shell  of  infinite  length. 
Consider  Figure  1-3 


,r 


Figure  1-3 

The  source  function  due  to  the  generators  of  the  cylinder  is  given  by  (1.4) 
and  dN/di  replaced  by  p^dOdr  and  R  by!  .  Therefore 


dp 


'Q\  -p/t^  (-ErRcose +r^  +  R^) 


(1.5) 


The  density  due  to  an  initially  uniform  expanding  cylinder  of  radius  a  is 
therefore 

_  ft  i»2  ^  -i  •>  r  ft  \ 

(1.6) 


p  =  2p, 


We  need  not  concern  ourselves  with  the  difficult  integrations  required  in 

(1.6)  provided  we  are  at  large  distances  (R;!>a)  from  the  cylinder,  for  then 

2 

i  =  R  and  (1.6)  reduces  to  (1.4)  where  dN/dl  s  . 
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Equation  (1. 6)  may,  by  a  simple  transformation  be  reduced  to 

2.2 

p=p^e'P*^'‘  /  dze'*  J^(i2Vy^  (I.6a) 

where  y  =  R^p  /t^. 

3  4 

The  integral  in  (1. 6a)  has  been  tabulated.  * 

3 

A.  D.  Wheelon  and  J.  T.  Robacher,  A  Table  of  Integrals  Involving  Bessel 
Functions,  The  R»W  Corp. ,  1954,  p.  57  Integral  No.  3.  108. 

4 

Math  Tables  and  Aids  to  Comp.  6,  40,  1952. 
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APPENDIX  n 


THE  INTERACTION  OF  AN  EXPANDING  GAS  CLOUD  WITH  A  SURFACE 
1.  Specular  Reflection 

a.  A  Source  of  Distance,  R,  from  an  Infinite  Plane 

Let  us  assume  we  have  a  point  source,  s,  a  distance,  R,  from  an 
infinite  plane.  At  t  =  0,  the  gas  molecules  are  emitted  and  travel  radially. 
Those  particles  that  strike  the  plane  are  assumed  to  be  reflected  specularly. 

We  determine  the  source  fvinction  for  this  problem.  Consider  Fig.  I1>1 


Figure  II- 1 

It  is  obvious  that  half  of  the  particles  emitted  by  s  are  reflected  from  the 
plane  and  upon  reflection  appear  to  emanate  from  s^^,  the  mirror  image. 
Because  of  this  symmetry,  the  distance  traveled  by  a  particle  from  s  to  the 
plane  to  a  point,  P,  is  exactly  the  same  in  magnitude  as  the  length  of  a 
ray  from  fo  P* 

Now  the  density  observed  at  a  point,  P,  on  the  source  side  is  com¬ 
posed  of  two  parts,  one  due  to  the  original  source  with  a  source  strength  dN, 
the  other  due  to  the  mirror  image,  s^^,  with  the  same  source  strength.  Thus 
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) 


dp 


(II.  1) 


+  dN 


Now  i  j  and  may  be  written  in  terms  of  R,  0  and  i  ,  as  in  Fig.  II- 1, 
but  we  shall  not  do  this  here.  Equation  (1)  may  be  used  as  a  source  function  to 
describe  the  evolution  of  any  gas  cloud  of  any  initial  configuration  in  front  of  a 
plane. 

The  method  of  images  employed  here  may  be  extended  to  more  compli¬ 
cated  cases  of  a  source  expanding  in  front  of  several  intersecting  planes.  But  in 
such  cases,  as  in  the  equivalent  electrostatic  case,  multiple  images  must  be 
considered. 


b.  A  Finite  Source  at  the  Center  of  a  Sphere 

A  source  is  released  at  t  =  0  at  the  center  of  an  evacuated  sphere  of 
radius  a.  The  particles  travel  radially  until  they  strike  the  wall  where  they 
are  reflected  and  then  return  to  the  center,  traveling  a  distance  2  a  all  together. 
The  particles  continue  to  travel  back  and  forth  from  the  center  to  the  boundary. 
Our  task  is  to  predict  the  evolution  of  the  gas  density  in  time  and  space. 

Let  us  assume  that  our  source  is  a  finite  gas  sphere  of  uniform  density 
initicilly.  It  is  at  the  center  of  a  much  larger  sphere,  so  much  larger  that  for 
most  of  the  region  within  the  sphere  the  evolution  of  the  gas  cloud  may  be 
described  as  though  it  came  from  a  point  source  containing  N  particles  [cf. 

Eq.  (9)]. 


Figure  II-2 
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We  now  ask,  what  particles  will  be  found  a  distance,  r,  from  the  center 
between  the  two  concentric  spheres,  r  and  r  +  d  r,  at  time,  t  (see  Fig.  II-2)? 
The  answer  is:  those  particles  with  velocities  Vj  =  r/t,  =  (2  a-  r)/t, 

Vj  =  (2  a  +  r)/t,  =  (4  a  -  r)/t,  =  (4  a  +  r)/t,  etc-  Therefore  we  may  write  the 

density  in  space  for  all  r  much  greater  than  the  initial  radius  of  the  source 
sphere 

P  ^  (2a-r)^  p/t^  (2  a  -  r)^ 

r 

2.  Diffuse  Reflection:  An  infinitesimal  source  interacting  with  a  diffusely 
reflecting  surface.  Let  us  consider  a  source,  s,  of  dN  particles,  released 
at  time,  t  =  0.  This  source  is  some  distance  in  front  of  a  convex  surface  which 
acts  as  a  diffuse  reflector  (see  Fig.  11-3). 


p(r»t)  = 


3/2 


1 

7 


Each  element  of  area,  d  A,  of  this  surface  is  assumed  to  reflect  any 
incident  particle  into  a  random  direction  within  a  solid  angle  2  it.  The 
reflected  particles  are  also  to  have  a  Maxwellian  velocity  distribution  with 
an  associated  temperature,  not  necessarily  the  temperature  of  the  surface. 
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We  are  to  determine  the  density  of  particles  everywhere  in  space  and 
time. 

Now,  the  number  of  particles  falling  on  d  A  per  unit  time  and  therefore 
emitted  per  unit  time  from  dA  is  just: 

N(t)  =  Y  dAcos  edp(t,R)  (11.3) 


where  dp  (t,  R)  is  described  by  £q.  (3).  We  now  proceed  as  in  the  development 
of  Eq.  (3). 

During  a  time  interval  t*  to  t'  4  At',  N(t')At'  particles  are  emitted  from 
dA.  Of  these  particles,  the  fraction,  f,  have  velocities  ranging  between  v  and 
V  -f  Av.  Now,  by  hypothesis  f  has  the  following  form: 


f 


2 

V 


Av  ♦ 


(11.4) 


Thus  at  a  time,  t  (greater  than  t'),  after  the  release  of  the  source,  particles 
of  number,  f  N  (t')At',  will  be  found  within  a  hemispherical  shell  of  radius,  r, 
and  thickness,  Ar.  It  is  obvious  that  r  =  v  (t-f)  and  r  +Ar  =  (v+Av)  (t-t'-Af). 
The  density  of  these  particles  at  r  is  therefore: 


dp 


Av 


2N(t')At' 

(^) 

Av 

(t-t')^ 

V  At '  +  Av  (t  -  t ' ) 

(II.  5) 


where  v  =  r/t-t'  and  dp^^  is  the  density  of  particles  reflected  from  the 
surface  with  velocity,  v. 

Now,  the  density  of  all  reflected  particles  at  r  and  at  time,  t,  is  obtained 
by  counting  particles  with  all  velocities.  The  velocities  for  particles  emanating 
from  dA  and  found  at  r,t,  range  from  a  minimum  of  r/t  to  infinity,  or  equivalently, 
t'  ranging  from  0  to  t.  The  density  of  reflected  particles  may  therefore  be 
obtained  by 
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(II.  6) 


Substituting  (II.  3)  for  N(t')  and  £q.  (3)  for  dp  we  obtain  the  formidable 


integral: 


3/2  .t 


j-  - y  Rd  A  cos  0  d  N  exp  ?  -  ^ 

t'  (t-f)  L  f 


(t-t-r 


(11.7) 


£q.  (II.  7)  can  obviously  be  generalized  to  the  case  of  a  finite  gas  source.  The 
density  at  r  due  to  the  whole  reflecting  surface  may  also  be  obtained  from  (11.  7) 
by  integrating  the  expression  over  the  surface.  The  density  of  all  particles, 
dp^  ,  will  be  given  by 

dp^  =  dp  +  dp^  (II. 8) 

where  dp  is  the  density  of  particles  due  to  the  free  expansion  of  the  source,  s. 

The  solution  of  even  the  simplest  problem  (e.g.  ,  a  point  source  in  front  of 
diffusely  reflecting  plane)  presents  formidable  analytic  difficulty  in  contrast 
to  the  case  of  specular  reflection.  No  such  solutions  are  presented  in  this 
paper. 

The  treatment  of  diffuse  reflection  from  a  concave  surface  requires  con¬ 
sideration  of  multiple  reflections  and  thus  increases  the  difficulty  of  an  already 
difficult  problem. 
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